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Silicene (a monolayer of silicon atoms) is a quantum spin-Hall insulator, which undergoes a topological 
phase transition into other insulators by applying external field such as electric field, photo-irradiation and 
antiferromagnetic order. We investigate the electronic and transport properties of silicene nanoribbons based on 
the Landauer formalism. We propose to determine topological phase transitions by measuring the density of 
states and conductance. The conductance is quantized and changes its value when the system transforms into 
different phases. We show that a silicene nanoribbon near the zero energy acts as a field-effect transistor. This 
transistor is robust though it makes use of the minimum quantized conductance since the zero-energy edge states 
are topologically protected. Our findings open a new way to future topological quantum devices. 



Silicene is a honeycomb structure of silicon atoms akin 
to graphene. Its experimental synthesis has opened a break- 
through in the study of silicene It has two salient features 
absent in graphene. One is the relatively large spin-orbit (SO) 
interaction, which enables quantum spin-Hall (QSH) effects 
to realize^. The other is its buckled structure, which enable us 
to apply different external fields between the A and B sub- 
lattices such as electric field and exchange field. As a re- 
sult we can externally tune the band gap of silicene and drive 
topological phase transitions from the QSH insulator to other 
insulators^"^. We have already shown how to make experi- 
mental observation of phase transition points with the use of 
diamagnetism^ and also optical absorption^. However, these 
methods may not be so practical. 

The most important graphene derivative is nanoribborU^— . 
The low-bias low-temperature conductance a of graphene 
nanoribbons has been shown to be quantized as 



cr = 4 (n + 1/2) (eV^), n = 0, 1, 2, • 



(1) 



both for zigzag and armchair edges. We interpret the formula 
to imply the electron-hole symmetry, the 4-fold degeneracy of 
each energy level associated with the spin and valley degrees 
of freedom, and the conductance quantum /h per chan- 
nel. It is interesting to examine the same problem in silicene 
nanoribbons because the spin-valley degeneracy is broken ac- 
cording to a specific pattern in each phase. 

Topological insulators are indexed ^^^^ by a set of two topo- 
logical numbers (C, Cg), where C and Cs are the Chem num- 
ber and the spin-Chem number modulo 2. There appear four 
types of insulators in silicene^'^'^^. They are the QSH, the 
quantum anomalous Hall (QAH), the spin-polarized quantum 
anomalous Hall (SQAH), and the trivial band insulators. The 
QSH effect is an analogue of the quantum Hall effect for spin 
currents instead of charge currents. The QAH effect is the 
quantum Hall effect without Landau levels, while the SQAH 
insulator has a hybrid character of the QSH and QAH insula- 
tors. The prominent feature of a topological insulator is the 
emergence of zero-energy edge states at half-filling which are 
topologically protected against perturbatioD^^. 

In this paper we calculate the density of states (DOS) and 
the conductance in zigzag silicene nanoribbons based on the 
Landauer formalism^^, and propose an experimentally bet- 
ter method to detect a topological phase transition by way 



of measuring them. Our first observation is that there are 
finite DOS due to the zero-energy edge states in a topolog- 
ical phase, while they are absent in the trivial phase. Con- 
sequently, the topological phase transition must be observed 
experimentally just by measuring the site-resolved DOS with 
controlled electric field, which can be achieved by spatially re- 
solved the scanning tunneling microscope/scanning tunneling 
spectroscopy (STM/STS). 

The particularly important quantity is the conductance due 
to the topologically protected zero-energy edge channels, 
which is the one experimentally observable. It is interesting 
that the helical edge of the QSH and the chiral edge of the 
QAH insulators have the same amount of conductance. Our 
result is summarized as 
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The conductance changes its quantized value when the system 
is transformed into different phases by tuning electric field. 
Because of this property a silicene nanoribbon may act as a 
field-effect transistor, where the conductance is quantized and 
topologically protected. 

Hamiltonian: The basic nature of silicene is described by 
the tight-binding Hamiltoniani2i2S, 



H 



{{1,3)) al3 
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where creates an electron with spin polarization a at site 
z in a honeycomb lattice, and (z, j) / ((^, j)) run over all the 
nearest/next-nearest-neighbor hopping sites. The first term 
represents the usual nearest-neighbor hopping with the trans- 
fer energy t = 1.6eV. The second term represents the effective 
SO interaction with Aso = 3.9meV, and Uij = +1 if the next- 
nearest-neighboring hopping is anticlockwise and Vij = — 1 
if it is clockwise with respect to the positive z axis. We have 
neglected the Rashba interactions since their effects are neg- 
ligibly small^ in general. Here, is the Pauli matrix for the 
the z component of the spin, whose eigenvalues are 5^ = ±1. 
We also use Sz =ti for indices. 
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FIG. 1: (Color online) Silicene nanoribbon in electric field, (a) It is 
decomposed into the device, right lead and left lead parts. The width 
is taken to be = 5. (b) Silicene consists of the A-sublattice and 
the 5-sublatice with layer separation 2^. The energy of the A sites 
(red disc) is lower than the one of the B sites (blue disc) in electric 
field Ez. The edge mode is localized along the A sites {B sites) 
of the up (down) outmost edge of a nonoribbon. See Fig|2]for the 
site-resolved DOS of the edge modes. 



There are two inequivalent Brillouin zone comers, called 
the K and K' points, as in graphene. The conduction and 
valence bands form two conically shaped valleys or cones at 
these points. We assign the valley index /i = ± to electrons 
at the K or K' point, and call them the Kr^ points as well. 
We also use r] = K^K' for indices. The low-energy effective 
Hamiltonian is given by the massive Dirac theory around the 
Kr^ point. The Hamiltonian ([3]) yields 



kyTy) + XSOV^ZCTZ 



(4) 



where = is the Fermi velocity with the lattice con- 

stant a = 3. 86 A. Here, is the Pauli matrix for the z compo- 
nent of the pseudospin representing the A and B sublattices, 
whose eigenvalues are = ±1. 

A great merit of silicene is that we can introduce various 
potential terms into the Hamiltonian by making advantages of 
its bucked structure. Eight commuting terms are possible in 
the Dirac Hamiltonian dU, 




FIG. 2: (Color online) Band structure, DOS and conductance of 
zigzag silicene nanoribbons for (a) the QSH insulator phase, (b) the 
metalic phase at the phase transition point, and (c) the trivial insula- 
tor phase. These phases are obtained by applying electric field Ez. 
The phase transition occurs at Ez = Ecr. The number of bands is 
2VK + 2 in the nanoribbon with width W. Here, the width is taken 
to be = 31, and only a part of bands are shown. The band gap 
is degenerate (nongenerate) with respect to the up (red) and down 
(blue) spins at =0 (Ez > 0). Van Hove singularities emerge 
in the DOS at the points where the band dispersion is flat. The site- 
resolved DOS of the up- spin state at the outmost A and B sites of 
a nanoribbon are shown by red curves in the insets. There are fi- 
nite DOS for the zero-energy edge states in the QSH insulator. The 
conductance is quantized by unit of /h. 



(5) 



where p, r = or 1. Each term has different symme- 
try properties. The coefficient of is the Dirac mass, to 
which four terms contribute. They are Hpqi. First, Hm is 
nothing but the SO coupling term with Am = Asq. Sec- 
ond, i^ooi is the staggered sublattice potential term^ with 
Aooi = ^Ez, where 2i is the separation between the A and 
B sublattices and Ez is the external electric field. Third, 
i^ioi is the Haldane term^\ where we set Aioi = A^^ by 
introducing photo-irradiation^'^^ with strength Xq. Finally, 
i^oii is the the staggered exchange tern>ii, where we set 
Aoii = AM = Ma — Mb by introducing the exchange fields 
Ma(b) to the A(5)-sublattice. 

We may write down the tight-binding term that yields the 
potential term Hpqi . The additional terms are^i^ 

- £Ez ti4^Cia + AM ^ ticl^(JzCia. (6) 



where t\ = ±.lfox i = A^ B. The Dirac Hamiltonian is 

Hrj = H^- lEzTz + vXnTz + AMa^r^. (7) 

The spin- valley dependent Dirac mass is given by 

A^^ = r]SzXso - £Ez + r]Xn + 5^ AM. (8) 

It may be positive, negative or zero. The silicene system is 
described by the Hamiltonian H + /S.H. 

Topological phases: Any insulating state is characterized 
by a set of two topological quantum numbers (C,Cs). Pro- 
vided the spin 5^ is a good quantum number, they are given 
by C = Cf + Cf + Cf + Cf and = i(Cf + Cf - 

— ), where C^^ is the summation of the Berry curva- 
ture in the momentum space over all occupied states of elec- 
trons with spin 5^ in the Dirac valley i^T^, and calculated^* 
as C^^ = |sgn(A^ ). All possible topological insulators 
are determined by the three parameters Ez, Xq and AM 
with respect to Aso- Possible sets of topological numbers 
are (0, 0), (2, 0), (0, (1, ^) up to the sign ±. They are the 
trivial, QAH, QSH, SQAH insulators, respectively. Note that 



3 




(0,1) £cr Ecr 



FIG. 3: (Color online) (a) Phase diagram in the {Ez,AM) plane 
and (b) conductance as a function of Ez at fixed values of AM. It 
takes a constant in one phase, and changes its value across the phase 
boundary. The unit is Aso for AM, and /h for the conductance. 
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FIG. 4: (Color online) (a) Phase diagram in the {Ez, Xn) plane and 
(b) conductance as a function of Ez at fixed values of Xq. It takes a 
constant in one phase, and changes its value across the phase bound- 
ary. The unit is Aso for Aq, and e^/h for the conductance. 



there are two-types of trivial band insulators, which are the 
charge-density-wave (CDW) type insulator and the antiferro- 
magnetic (AF) insulator 

A topological phase transition occurs when the band gap 
closes, or A^^ = 0. Let us review^ the topological phase 
transition along the Ez axis, where the Dirac mass is given 
by A^^ = r]SzXso — ^Ez- The condition A^^ = implies 
Ez = ±E^, with E^, = Xso/i. It follows that (C, Cs) = (0, 0) 
for \Ez\ < E^, and (0, \) for \Ez\ > E^,. We have given the 
band structures at Ez = 0, E^r and 2Ecr in Fig|2l It is the 
characteristic feature known as the bulk-edge correspondence 
that zero-energy edge modes emerge in topological insulators 
with (C, Cs) / (0, 0) and that these zero-energy edge modes 
are topologically protected against perturbation 

We may easily construct the phase diagrams by solving 
A^^ = 0. Those in the (Ez.AM) space and the {Ez.Xq) 
space are given in Fig sl3] and |4] together with the pair of topo- 
logical charges (C,Cs), respectively. All possible topological 
insulators are found in these phase diagrams. The band struc- 
tures at typical states are found in Figsl2] and E] where the 
breakdown of the spin- valley symmetry is manifest. We pro- 
ceed to characterize each phase by its characteristic pattern of 
the DOS and the conductance of a zigzag nanoribbon. 

DOS and conductance: The natural framework for trans- 
port calculations in nanoscopic devices is the Landauer 
formalism^^. We consider a zigzag silicene nanoribbon di- 
vided into three regions [FigJT]: the device region, the left 
lead and the right lead. The size of the device region is actu- 
ally irrelevant due to the ballistic transport property. 

In terms of single-particle Green's functions, the low-bias 
conductance cr{E) at the Fermi energy E is given by^ 

cj{E) = {eyh)Tr[ri^{E)GUE)T^{E)Gj,{E)], (9) 

where rR(L)(£;) = i[^R(L){E) - with the self- 

energies ^l{E) and Er(£^), and 

G^{E) = [E-Hj,- El(^) - ^r{E)]-\ (10) 

with the Hamiltonian i^D for the device region. The self- 
energy I;l(r)(£^) describes the effect of the electrode on the 
electronic structure of the device, whose real part results in a 
shift of the device levels whereas the imaginary part provides 
a life time. It is to be calculated numerically^—. 




FIG. 5: (Color online) Band structures of the SQAH and QAH insu- 
lators. Up- spin (down- spin) states are illustrated in red (blue). There 
are two channels in the QAH insulator but only one channel in the 
SQAH insulator that contribute to the conductance at half-filling. 

The total density of states (DOS) reads 

p{E) = -^-ilmTrGD(^), (11) 
while the partial density of states at i site reads 

p,{E) = -n-'lm[Gj,{E)u], (12) 

in terms of the Green function Gd{E) of the device. 

We have calculated the DOS p{E) and the conductance 
cr{E) of a nanoribbon as functions of the Fermi energy E, 
which is controlled by doping. We give the results at electric 
field Ez = 0, Ecr and 2Ecr in Fig 121 A van Hove singularity 
occurs in the DOS at the point where the band dispersion is 
flat. As E increases beyond the point, the Fermi level crosses 
a new band. A new channel opens and contributes to the con- 
ductance by e'^/h for each spin and valley. It is clearly ob- 
served that the edge channel connects the tips of the Dirac 
cones with the same spin at the K and points. 

We have also plotted the site-resolved DOS Pi{E) of the 
up-spin states at the outmost A and B sites of a nanoribbon 
by red curves in the insets [Fig 121. They represent degenerate 
zero-energy states at Ez = 0. As we have explained in Fig[T] 
the energy of the A and B sites become different for Ez ^ 0. 
It results in the downward (upward) shift of Pa{b){E) along 
the edge as Ez increases. They are separated completely, and 
zero-energy states disappear for Ez > Ecv 

The zero-energy edge channel of a topological insulator is 
particularly important at half-filling, because it is topologi- 
cally protected. We have calculated the conductance at half- 
filling by increasing the external field Ez in the {Ez, AM) 
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space and the (Ez^Xn) space. First, we increase from 
Ez = Sit AM = 0, |Aso, f Aso in the (E^.AM) space 
[Fig El- Initially the conductance reads a = 2e^//i for 
AM = 0, where the system is in the QSH phase. It reads 
a = e'^/h, when the system is in the SQAH phase. Its band 
structure is shown in FigI3a), where the zero-energy edge 
states account for the conductance a = e'^/h. The conduc- 
tion becomes zero as Ez increases and the system becomes 
the trivial AF insulator. 

We confirm these observations by investigating the same 
problem in the {Ez^Xn) space [FigJU. The system is in the 
QAH phase for > Aso at Ez = 0, where a = 2e^ /h. We 
illustrate the band structure of the QSH insulator in FigJSjb), 
where the zero-energy edge states account for the conductance 
cr = 2e^ /h. The edge channel is helical (chiral) in the QSH 
(QAH) phase, but both of them transport the same amount of 
electric charges when the current is fed. We summarize the 
conductance in each topological insulator as in (|2]). 

Field-effect topological quantum transistor: The con- 
ductance is quantized in silicene nanoribbons. The simplest 
system is provided by applying electric field only^, where 
quantized conductance changes from 2 to at the critical elec- 
tric field £^cr. This means the system acts as a transistor where 
"on" state can be switched off to "off" state by applying elec- 
tric field. This transistor is "quantum" since the conductance 
is quantized, which is highly contrasted with the ordinal tran- 
sistor, where the conductance is not quantized. Furthermore 
the conductance is topologically protected because the zero- 
energy edge state is topologically protected. Namely the con- 
ductance is robust against impurities due to its topological sta- 
bility. Consequently we may call it a field-effect topologi- 
cal quantum transistor. This is the most energy-saving device 
since it utilizes the minimum conductance. 

We are able to design a three-digit quantum transistor by 



attaching antiferromagnet. Namely the conductance changes 
in three steps 2, 1, with increasing electric field when 
I AM I < Aso- When |AM| > Aso, the conductance changes 
in three step 0, 1,0. It acts as a three-step transistor, where the 
system is first in the "off" state, then become "on" state and 
finally become "off" state with the increase of electric field. 

Conclusions: We have analyzed the DOS and the conduc- 
tance in silicene nanoribbons. There are finite DOS due to the 
zero-energy edge states in the topological phase, while they 
disappear in the trivial phase [FiglJl. Local DOS measure- 
ment is a direct evidence of the existence of the edge states, 
which can be achieved by spatially resolved STM/STS. This 
must be the most efficient way to make an experimental ob- 
servation of a topological phase transition. Furthermore we 
can determine the band gap by measuring the DOS. A precise 
measurement is possible owing to the van-Hove singularities 
present at the tips of the conduction and valence bands [Fig El • 

Conductance measurement is also a direct method to ob- 
serve a topological phase transition. We have proposed a field- 
effect topological quantum transistor with the use of the zero- 
energy edge state of a silicene nanoribbon. This could be a 
basic component of future topological quantum devices. 

In passing we address the problem how narrow the nanorib- 
bon can be. The penetration depth of the zero-energy edge 
state has already been shown to be as short as the atomic 
scale^^ in zigzag nanoribbons. Hence we may use a quite 
narrow nanoribbon to detect and make use of the topological 
properties of silicene. 
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